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Abstract 

The critical buckling characteristics of hydrostatically pressurized complete spher- 
ical shells filled with an elastic medium are demonstrated. A model based on small 
deflection thin shell theory, the equations of which are solved in conjunction with 
variational principles, is presented. In the exact formulation, axisymmetric and 
inextensional assumptions are not used initially and the elastic medium is mod- 
elled as a Winkler foundation, i.e. using uncoupled radial springs with a constant 
foundation modulus that is independent of wave numbers of shell buckling modes. 
Simplified approximations based on a Rayleigh-Ritz approach are also introduced 
for critical buckling pressure and mode number with a considerable degree of ac- 
curacy. Characteristic modal shapes are demonstrated for a wide range of material 
and geometric parameters. A phase diagram is established to obtain the requisite 
thickness to radius, and stiffness ratios for a desired mode profile. The present 
exact formulation can be readily extended to apply to more general cases of non- 
axisymmetric buckling problems. 
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1. Introduction 



The analytical study of the structural behaviour of spherical shells is of great 
importance in the fields of not only civil, mechanical and aeronautical engineering 
but also nano science and biomechanics; Notable examples include pressure ves- 
sels, spherical honeycombs (\], carbon onions and spherical viruses ^ and 
so on. In the engineering research fields, some pioneering theoretical works on 
the elastic instability issues of empty complete spherical shells were conducted in 
1960s by Thompson ||3], Hutchinson [S] and Koiter [Q]. 

In recent years, analytical works which considered an interaction between a 
spherical shell and an internal elastic medium have also been conducted. The 
interaction effect leads to novel buckling patterns which depend on the stiffness 
and thickness to radius ratios. Recently for example, Yin et al. [lZO successfully 
demonstrated the stress-driven buckling patterns in spheroidal core/shell struc- 
tures, where the core implies an elastic medium, by using the finite element ap- 
proach. Interestingly, the authors insisted in that work [7 ] that shapes of many nat- 
ural fruits and vegetables can be reproduced by anisotoropic stress-driven buckles 
on the spheroidal core/shell system. Moreover, the morphology of a pressurized 
spherical shell can give clues to the structure-property relationship of spherical- 



9|1, core/shell semiconductor 
The last system refers to a 



shaped nanostructures such as carbon onions m 
nanoparticles lIlOll . and the "nano-matryoshka" [IH 
concentric multilayered structure comprising various metallic and dielectric ma- 
terials, showing great tunability in the optical response; it is thus of interest to 
explore how pressure-induced deformation affects the electromagnetic properties 
of the nanoshells and their resonant frequencies. 

Purely from a structural mechanics perspective, Timoshenko |fl2!], and Fliigge 
[|l3n introduced the formulation for buckling of a complete spherical shell in their 
respective books that are now regarded as classics in the field. In these books, 
the exact approaches used to solve the hydrostatically pressurized buckling of 
a complete spherical shell without an internal elastic medium were introduced. 
However, both of these formulations were based on the axisymmetric assumption 
and formulations for more general cases including non-axis ymmetic deformations 
were not included. More recently, Fok and AUwrite [|l4ll analysed the elastic 
axisymmetric buckling behaviour of a complete spherical shell embedded in an 
elastic material and loaded by a far-field hydrostatic pressure. In that study, the 
energy method in conjunction with a Rayleigh-Ritz trial function was used for 
simplicity but the validation of the obtained results was not discussed in detail with 
buckling deformation modes being omitted. Hence, to our knowledge, no general 
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non-axisymmetric formulation, in conjunction with exact methods of solution, has 
been developed for the buckling behaviour of a complete spherical shell with an 
intemal elastic medium thus far. 

Moreover, much attention has been given to the structural morphology of the 
core/shell structure. Some authors have recently demonstrated the cross-sectional 
morphology of carbon nanotubes embedded in an elastic medium lIlSL II6L 11711 . 
These results clearly show that interactions between shells and cores lead to some 
novel wavy- shaped buckling deformation patterns. 

From the background described above, the buckling properties of hydrostati- 
cally pressurized complete spherical shells filled with an elastic medium is demon- 
strated currently. The elastic medium is modelled as a Winkler foundation, i.e. 
with uncoupled radial springs and a constant foundation modulus. An exact ap- 
proach for solving the developed equations based on the formulation without us- 
ing the axisymmetric and inextentional assumptions is presented. This approach 
therefore avoids any discussion about the validity of the solution and allows to 
model to be extended to cover more generic non-axisymmetric cases with relative 
ease. The analytical results are presented using a phase diagram and illustra- 
tive buckling modes. In addition to this, simplified approximations based on a 
Rayleigh-Ritz approach are also introduced for the critical buckling pressure and 
the corresponding mode number. Comparative studies between the exact and sim- 
plified approaches are conducted to validate the approximation and it is shown that 
the approximate formulations enables sufficiently accurate values to be obtained. 



2. Complete spherical shell model 

The critical pitchfork bifurcation phenomenon of the hydro statically pressur- 
ized complete spherical shell is investigated, as shown in FigH] The spherical 
shell is constructed from an homogeneous and isotropic linear elastic material 
with Young's modulus E and Poisson's ratio v, which is filled with an elastic ma- 
terial that is modelled as a Winkler foundation, i.e. with uncoupled springs in the 
radial direction and a constant foundation modulus kf. For a complete spherical 
shell with radius a and thickness h, spherical angular coordinates in the latitude 
and meridian directions 9 and ^ are used. Displacement functions: u, v and w are 
in the 9, ^ and the (outward) radial directions, respectively. 
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Figure 1 : Hydrostatically pressurized spherical shell filled with an elastic medium. 



3. Formulation 

3.1. Exact approach 
3.1.1. Energy formulation 

The following analysis is based on classical small deformation theory of thin 
shells [Il8l.ll9ll. The total potential energy V is expressed by the sum of the strain 
energy and work done by external load as: 

V = UM + UB + UF + n, (1) 

in which L^m is the membrane (in-plane) strain energy term, where: 

Umg = 1 [ [ NeEea^ sin^ d^ dG, 
2 J(p Je 

UM0e= J J^N^^B^ea^^i^^d^dG, (2) 
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Ub is the bending (out-of-plane) component, where: 



Ub = UBth + UBO +Ub09, 



<j) .16 



M^X^a'-?,m(^ d0, 

\ ! ! MeXea^^m^ d(^ d0, 
2 Je 

i7/7 is the strain energy term due to a Winkler foundation, thus: 

Up = ]r ( f a^kfw^ sin d) dG, 
and is the potential energy of the applied pressure: 

Q.= I I a^pwsm^ d<^ dO. 



(3) 



(4) 



(5) 



Based on the assumptions of thin shell theory [|l8ll . the strain-displacement 
relations can be expressed as: 

1 



-(v,^+w), 



coi(pv + — — -+w 



sin^ 



1 / V, 



a \sin^ 



COt(^M + M,(^ 



1 



a 



J8, 



1 



e 



a \ sin^ 
1 



+ u 



X(l> = ■^[-w,(j„j> +v„ 



Xe = 



1 



— —+u,e+coi(l>[v-w,(^ 



sin 



sin(^M,(|, +v,0 —2w,0o +2cot^w,e —cos^u 
2a^ sin ^ 



(6) 
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where denotes differention with respect to x, and the constitutive relations 
are as follows: 



^0 =C(£^ + V£e), Ne=C{£e + v£^), 
1 -V 

1 -V 

where C = Eh/{\ — v^) is the membrane stiffness and D = Eh^ /\2{\ — v^) is the 
flexural rigidity of the spherical shell, respectively. Substituting the expressions 
in Eqs.® and (|7]) into Eqs.©-© gives the energy expressions in terms of the 
displacements u, v and w. 

3.1.2. Fundamental state 

When uniform hydrostatic pressure p is acting on the spherical shell, only the 
inward radial static displacement wq can occur in the pre-buckling state. Hence, 
the stress-strain relationship of the spherical shell is assumed to be linear up to the 
point of instability. In the fundamental pre-buckling state, the potential energies 
y(*^) are given by: 

y (0) = (0) + (0) + f/f + ^^(0) = f f p{0) d0 (8) 

Jtj) Je 

where 

= C(l + V) / [wld^ dO, C/f = 0, 
J(l> Je 

[/(o) = \ I f ah fw^ sin ^ d0 d0, 
2 J(j) Je 

^^(0) = / / a^pwosivKj) d<^) dd. (9) 
J(l) Je 

For an equilibrium state, the first variation of the total potential energy V must 
equal zero. This condition gives the following static displacement under uniform 
hydrostatic pressure: 

2Eh + a^kf{l — v) 
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3.1.3. Critical buckling analysis 

To obtain an expression for the second variation of the total potential energy, 
the following infinitesimally small increments are defined: 



M = Ml, V = Vl, w = wo + wi, (11) 

which correspond to buckling displacement modes. The expression for the po- 
tential energy due to the external force is a linear functional of the displacement 
components and makes no contribution to the second variation expression, which 
is — [|l8n . Consequently, the second variation of the total potential energy 
becomes 

5^U = 5^UM + 5^UB + 5^Uf= [ [ Fd(j)de, (12) 

Jd J 9 



where 



5^Um —S^Ustretch + ^^Usheari 

-2tt _Cwo(1 + v 



S^Ustretch — [ [ [(vi+wi,w,)^ + (mi -Wi,0CSC0)^]sin0 d^ d0 

a J^Je 

+ C [(vi + -wi,ecsc^)^] d^ d0, 

Je 

5^Ushear= ^ / / (cot + CSC ^vi ,e +Mi )^ sin d0 d0, (13) 



JO 



5^UB=5^Ube„d + 5^Utwist, 
5^f/w=^ / / {(cOt^Mi +CSC^Vi,e-CSC^0Wi,00-COt0Wi,^)^ 

a^J(j)Jd 

(vi„^-|-wi„^,^)[vi„|,+wi,^^+2v(cot0vi +csc0Mi,e-csc^^wi,00 
— cot^wi, tj) )]}sin^ dO, 

S^Utwist = ^^\^2^^ ^^[2(1 + v)wi(viCOt0+Wi +Mi,e+Mi,0) 
+ (2VVlCOt^+2VMl,0+Vl,0) 

+ (vicot^ + Mi,0csc^)^] sin^ d^ d0, 

(14) 

d^UF= I [ a^kfw\sm^d(l)de. (15) 



According to the Trefftz criterion, the buckling equations can be obtained by 
introducing F in Eq.(fT2l) into the Euler-Lagrange equations with the calculus of 
variations. The Euler-Lagrange equations in this case are as follows fi 1 811 : 



dF d f dF \ d ( dF 



du\ dcj) \dui_(j) J do \dui^e 
dF d f dF \ d f dF 



0, 



dv\ d({) \dvi_(j) J do \dvi 
dF d f dF \ d f dF 



dw\ d^ \dwi ^ J do \dwi q 

^4f^V^f^V4r^Ho, (16) 



(9(^2 \dwi^(j,^ J d^dO \dwi^^Q J dO^ \dwi^QQ 

Substitution of a solution of the form 

Ml = ^m(^) sinn0, 
vi = ^v(0)cosn0, 

wi = ^vi>(0)cosn0, (17) 

where n is a positive integer, into Eqs. (fT2)) -(fT5l) and some rearrangement gives the 
following ordinary differential equations for an arbitrary n 

„»|i .^(1 + v)(«d„^, = 0, (18) 

d(p^ dip' dip' a 

(i)d'u ,(i)d'v (i)d'w WO/, X . ,/ dw, „ 

wq , X . , / <9v , nu d^w dw , \ ^ 

+ — (1 + v) sm0 — +vcot0 + ^— - ^-y - ^cot(^ H = 0, 

(20) 



in which the implied summation over i cover the range of non-zero coefficients, 
which are described in the Appendix. These governing differential equations 
(fT8l)-(|20l) are difficult to solve exactly by only assuming periodic functions in w, v 
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and w directly. However, the definition of the following special linear differential 
operator enables the solution of these equations with an exact approach: 

//(...) = (...)"+ (...y cot +(...) (2-^) , (21) 

V sin ^ / 

where primes denotes differentiation with respect to . From the definition of H 
in Eq.(l2TI). the following relations are readily obtained: 



HH{...) ={...)"" + 2{...y"cot(^ + {...)" (^3 

l+2n2~ 



+ (...)'cot0 4 + 



sin 



sin 

Sin f sin f 



H'{...) = (...)'"+ (...)"cot0 + (...)'cot<^, (2 - 1±^) 

V sin / 

In addition, a new variable m is defined with regard to u by 



+ 2„^(...)^^. (23) 



sin 



By making use of the relations of Eqs.(l2TI)-(l23l) and the new definition of Eq.(l24l). 
Eqs.(fT8l) and (fT9l) can be rewritten respectively as: 



u = u . 



(24) 



d (\-v 



H{u) 



1+v vs'm^ +nu 

— ^ — 7 

2 sin (j) 



^.^ . vsin^+nu 
— 2(1+ k)n — cos (j) 



sin 



-{l+k){l+v)nw + knH{w)- {u + nw) = 0, 



(25) 



■\+k) 



^-^^{vsin^ + nU)" -\ ^//(vsin(^) — ^ ^ ^ cot0(vsin0 



^1 +v)w'sm^ 



■k%m^H\w) _ ^o(^ + sin<^(v-w ) = 0. (26) 



Differentiation of Eq.(|26l) with respect to ^ and then multiplying it by sin(^/R 
leads to the following: 



sin(^ f 1 + V . 



+ ■ 



3-v 



sin 



2 ,vsin(^+ni2)'" + ^— //'(vsin(^) 



(vsin<^ — cot^(vsin(^ 



+ (1 + v)(w"sin<|) +w'cos0) I -/t^^[cos(|)//'(w) + sin0//"(w)] = 0. (27) 

J " 
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By adding Eq.(|271) to Eq.(|25l). some rearrangement gives: 



( 1 + fc) [//(^) - ( 1 + v) - //(w) + 2w)] 

-H{w)+2w)=0. 

(28) 



-k {HH{w) - 2//(w) ) - ^ ( 1 + v) - // (w) + 2w) = 0. 



in which: 

^^(vsin0+n^ (29) 
sin0 



Similarly, Eq.(l2()l) can be written using the operator H and the new variable ^ 
defined by Eq.(|29l) as follows: 



( 1 + yt) ( 1 + v) + 2w) + k[HH{w) - (3 + v)H{w) -H{^>)\ 

+ ^w + (l + v)— (^-//(w)+2w) =0. (30) 
C a 

Equations (|28] ) and (|30l ) are the governing equations to be solved. Now, the 
solutions of m and w of Eq.(|28]) and (|30l ) is assumed to be in the form 

oo oo 

^= A„P^(cos(^), w= 5„,i',„(cos(^), (31) 

m=0 )7j=0 

where A„, and 5„, are the constant deformation amplitudes, and the spherical har- 
monic Pm(cos^) is a series of Legendre functions of degree n which satisfies the 
following Legendre differential equation, 

P,;; + P'„, cot + m(m + 1 )P,„ = 0. (32) 

By using the operator// defined in Eq.(l2T]). the following relation can be obtained: 

H{P,n) = -X„rP,n, (33) 

where 

A„, = m(m+l)-2, (34) 

and 

HH{P,n) = -XlP,n. (35) 
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Substituting the series of Eq.(|3T|) into Eqs.(|28]) and (l30l ) and using Eqs.(l33l) and 
( 1351 ). we obtain: 



(c2^U„, + c),!^''5„,)/'„,(cos(^) = 0, 



.(12 



m=0 



^(cL 'A„^ + clJ 'B„,)P,n{cOS^) =0, 
m=0 



(36) 



in which 



(11) 



J12) 



.(21) 



.(22) 



Wo 



;i+^)(A,„ + l + v) + — (1 + v 



Wo, 



=(1 +/:)(! + v)(A,„ + 2)+fc(A,f, + 2a,„) + —(l+v)(l+A,„), 



=(i+fc)(i + v)+a„,+— (i+v), 



a^kf Wo 



=2(l+fc)(l + v)+a,„(A,„ + 3 + v) + — ^ + ^(l + v)(A,„ + 2). (37) 

C a 



Hence, homogeneous linear equations for A,„ and 5„, are obtained for all m as 
follows: 
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12 



^(21) ,.(22) 



(38) 



This is a standard eigenvalue problem for which nontrivial values are obtained 
when the coefficient matrix in Eq.(l38T) becomes singular. The corresponding pres- 
sure p is thus calculated from the determinantal equations: 



det 



^(12) 



^(21) ,.(22) 



0. 



(39) 



By solving Eq.(|39l). the pressure p with the corresponding mode m for several 
combinations of geometric and material constants can be obtained as an eigen- 
value. For the particular values of the constants, the minimum value of the eigen- 
values is the critical buckling pressure Pcr- It should be noted that the above equa- 
tion is independent of the mode number n. This shows that only axisymmetric 
modes can occur in this problem despite the inclusion of the Winkler foundation 
term. 
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3.2. Simplified approach 

Next, we consider the simplified formulation by the Rayleigh-Ritz approach. 
Here an inextensional and axisymmetric buckling deformation of the shell are 
assumed with shear and torsional strains being neglected for simplicity. Thus, u, 
e^Q, and XipQ are assumed to be zero. 

The buckling deformation forms are given as follows: 



Vl = \/l ^ V,n 

m=2 



dP,n{c) 

dc 



w 



1 = 52 "^mPmic) 



(40) 
(41) 



m=2 



where c = cos0. 

The first variation 5Um of the membrane (in-plane) strain energies is given by 
5Um = 1tzC{1 + v) J Wo (v'l +vicot^+2wi)sin0(i0. (42) 



The integrand in Eq. (|42l) . 



v[ +Vi cot^ + 2wi = 0, 



(43) 



gives the inextensional buckling condition and when the assumed deformation 
functions in Eq.(l40l) and (|4TI) satisfy the condition, the following relation can be 
obtained 

m(m+ 1) _ 



The second variation of the strain energies are 

5^V = 5^Um + 5^Ub + 5^Uf = 2k I Fd(j), 

J(h 



(44) 



(45) 



in which the each term can be expressed by the strain-displacement relations and 
Eqs.dini) and dH]) as: 



2(m- l)m(m+ l)(m + 2)(l - v) 



m=2 



2m +1 

Wo (m- l)2m(m+ l)(m + 2)2(l + v) 
a 2m + 1 



~2 

V™ ■ 



(46) 
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d^UB 



nD ^ (m-l)2m(m+l)(m + 2)2[m(m+l)-l + v]_2 



m=2 



The stability criterion 



m=2 
2 



2m +1 

s:2tt 2, m2(m+ 1)2^2 



dv,n 



0, 



(47) 
(48) 

(49) 



gives the following critical pressure 

2Eh + kfa^{l-v) j 2Eh + {m- \){m + 2) 
^" ~fl£/z(m-l)2(m + 2)2 j r+v 

D(m- l)2(m + 2)2[m(m+ 1) - 1 + v] 2 



(50) 



+ 



+ a kfm{m+ 1) >. 



Equation (1501) has a minimum value when 

dpcr 



(51) 



In such spherical shell buckling, m is usually a large integer. This allows us to 
approximate (m — l)2(m + 2)^ ^ m'^, m(m + 1) ~ m^, and so on. By applying the 
approximations to Eg. (1501) and substituting the resulting equation into (|5TI) . the 
mode number m^r associated with the lowest critical pressure and the correspond- 
ing critical pressure per can be obtained as follows: 



2Eh + ka^{l-v) 
aEhmi^ 



la^[2Eh + a^k{\+v)] 
D{\ + v) ' 

2Eh + ka^(l + v) DnP- 



1+v 



+ ■ 



~ 2 
111 ■ 

cr\ cr 



1 + v) 



(52) 



(53) 



4. Analytical results and discussion 

4.1. Comparison of buckling pressures between exact and simplified approaches 

Figure |2] shows the variation of the nondimensionalized critical buckling pres- 
sure Per/ Po, in which po is the critical pressure for the empty complete spherical 
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Figure 2: Comparison of pc- {a/h= 100). 
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Figure 3: Relative error of pcr {ajh = 100). 
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Figure 4: Comparison of nicr {a/li = 100). 
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Figure 5: Relative error of nicr {ajh— 100). 
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shell, which is given by: 



2E fh^ ^ 



"-75?!=^'^' ■ ''''' 



The solid lines are the exact values obtained from Eq.(|39i) and the dotted lines 
correspond to Eq.(|53]l from the Rayleigh-Ritz analysis. As can be easily seen in 
this figure, the value of p^r/ Po increases with increasing stiffness ratio akf/E and 
the exact and simplified curve agree with each other. Fig. |3] gives the relative 
error between exact and simplified values, which is defined by 

/'cr(simplified) /'cr(exact) /ec\ 
Pcr{exacl) 

The relative error converges with the increase of akf/E. 

The comparison of the critical buckling mode number nicr obtained from the 
exact and simplified approaches is shown in Fig.|4] Larger values of akf/E cause 
higher buckling modes, as in the case of beams on an elastic Winkler foudation 



[|20|] . The relative error 

'^cr(simplified) '^cr(exact) 



Mr 



(56) 

'cr(exact) 

is given in Fig. [5l We can also see the convergence behaviour with the increase of 
akf/E. These comparisons in Figs. |3] and |5] confirm the validity of the proposed 
simplified formulations for pcr and m^T especially for higher stiffness ratios. 

4.2. Eigenmodes 

Figure |6] shows a phase diagram of the buckling mode number m. This figure 
gives significant information on the various axisymmetric buckling modes de- 
pending on the values of akf/E and and a/h. It is observed from this figure that 
larger akf/E and a/h values favour buckling modes with higher mode numbers. 
We can also find that when the spherical shells are relatively thin with larger a/h 
values, the effect of the foundation modulus becomes more significant when com- 
pared with the case for the relatively thick shells with smaller a/h values. It should 
be noted that The obtained mode number m=18 for the condition of a/h — 100 and 
akf/E = 0. 1 is same as that found by Koiter |@]. 

Figure |7] illustrates the buckling eigenmodes for the empty complete spher- 
ical shell of a/h = 50. The characteristic wavy-shaped axisymmetric buckling 
deformation with the mode number m = 12 can be found in this case. 
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a 



h. 

Figure 6: Phase diagram of the buckling mode number m. Various axisymmetric buckling modes 
can be found depending on the values of akj/E and and a/h. 

On the other hand, Figs.[8]and|9]show the eigenmodes for the "filled" spherical 
shell {akf/E = 0.1) with a/h=50 and 100, respectively; the axisymmetric buck- 
ling deformation of the mode numbers m = 20 and 33 are seen in these figures. As 
can be understood in the developed formulation, non-axisymmetric modes cannot 
occur under the condition of constant foundation modulus. However, the cross- 
sectional shape for Fig. |9]is moderately different from that for Figs. |7]and[8l 

The results shown in Figs. suggest that the eigenmodes change corre- 
sponding to the mode number m. When m is even, the corresponding deformation 
of the spherical shell is symmetric about the equator, as shown in Figs. |7] and [8l 
On the other hand, when m is odd, the deformation is "anti-symmetric" about the 
equator, that is, an inward deformation at the north pole is accompanied by an 
outward deformation at the south pole. This is due to the property of the geom- 
etry of the sphere and the Legendre functions, not whether the shell is empty or 
otherwise. 
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Figure 7: Buckling modes of a hydrostatical pressurized "empty" complete spherical shell with 
a/h — 50 and akf/E ~ 0. (Left) 3-D view. (Right) Cross-sectional view through a great circle. 
The buckling mode number m— 12 is found. 




Figure 8: Buckling modes of a hydrostatical pressurized "filled" complete spherical shell with 
a/h — 50 and akf/E = 0.1. (Left) 3-D view. (Right) Cross-sectional view through a great circle. 
The buckling mode number m = 20 is found. 




Figure 9; Buckling modes of a hydrostatical pressurized "filled" complete spherical shell with 
a/h = 100 and akf/E = 0.1. (Left) 3-D view. (Right) Cross-sectional view through a great circle. 
The buckling mode number m = 33 is found. 



18 



5. Conclusions 



The characteristic buckling modes in hydrostatically pressurized complete spher- 
ical shells filled with an elastic medium have been demonstrated. A theoretical 
formulation based on small-displacement thin shell theory has derived governing 
equations of equilibrium that have been solved using an exact methodology. In 
addition, the simplified formulations for estimating the critical pressure and the 
mode number are proposed with a Rayleigh-Ritz approach. The formulations ob- 
tained from the Rayleigh-Ritz methodology is shown to give sufficiently accurate 
results compared with the exact values. 

From the formulation developed currently only axisymmetric eigenmodes can 
occur in spite of adding the Winkler foundation term; critical modes that are sym- 
metric or anti-symmetric about the equator may be determined depending on the 
combination of the foundation stiffness and the radius to wall-thickness ratios. 
A phase diagram has been established to obtain the requisite values of a/h and 
akf/E for observing a corresponding buckling mode. 

This work is a fundamental investigation that has been developed with future 
studies in mind. These will address the post-buckling behaviour [12 iL 12211 and 
non-axisymmetic modes which may occur when a foundation modulus is not a 
constant value and a function of the wave number of the shell deformation. The 
Rayleigh-Ritz approach developed currently is more likely to be employed for the 
nonlinear studies given that it allows for a simpler solution methodology; it has 
been shown in the present work that it provides an excellent approximation to the 
linear case. This will also provide an insight into the electrical properties of carbon 
onions, i.e., quasi- spherical nanoparticles consisting of concentric graphitic shells. 
They possess a high specific surface area (~ 500 m~ per gram) with no porosity 
and high field-emission stability owing to the endurance of carbon-carbon atomic 
bonds. These two features are appealing in the design of electrochemical capac- 
itors n23n and point electron sources [|24ll . in which the radial buckling of outer 
layers may lead positive effects in operation. It was also suggested in ||8|] that the 
electronic structure of carbon onions is sensitive to the radii of concentric layers 
and their separations. Hence, their asymmetric buckling under pressure should 
cause a significant change in the electronic states, where the mechanical deforma- 
tion of the underlying carbon structure induces a effective electrostatic potential 
that exerts on the mobile electrons on the layers [1251 \26n . The present findings 
will be helpful for addressing the deformation-driven carbon-onion behaviours 
from both academic and practical viewpoints. 
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Appendix A. The coefficients of the governing equations 

The coefficients of the goveming equations (fT8l)-(l20l) are written as follows: 
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